{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "from utils import *\n",
    "from mpi4py import MPI\n",
    "\n",
    "comm = MPI.COMM_WORLD\n",
    "rank = comm.Get_rank()\n",
    "size = comm.Get_size()\n",
    "\n",
    "maxTimePerSave = 5e-2\n",
    "\n",
    "save_dir = \"data/DME_compare/\"\n",
    "mech_arr = [{\n",
    "                'name':\"DMEzhao\",\n",
    "                'path':\"mech/DME/DMEzhao/chem.xml\"\n",
    "            },{\n",
    "                'name':\"DMEsk42\",\n",
    "                'path':\"mech/DME/DME_drg/DMEsk42.cti\"\n",
    "            },{\n",
    "                'name':\"DMEsk43\",\n",
    "                'path':\"mech/DME/DME_drg/DMEsk43.cti\"\n",
    "            },{\n",
    "                'name':\"DMEsk47\",\n",
    "                'path':\"mech/DME/DME_drg/DMEsk47.cti\"\n",
    "            },{\n",
    "                'name':\"DMEsk30\",\n",
    "                'path':\"mech/DME/DME_drg/DMEsk30.cti\"\n",
    "            }]\n",
    "T_arr = [700., 800., 1000., 1200., 1600.]\n",
    "P_arr = [1.*ct.one_atm, 10.*ct.one_atm, 20.*ct.one_atm]\n",
    "phi_arr=[0.5, 1.0, 1.5]\n",
    "\n",
    "props = {'T': 1000.,                # Temperature\n",
    "         'P': 30.*ct.one_atm,       # Pressure\n",
    "         'phi': 1.0,                # equivalence ratio\n",
    "         'air':'O2:1.0,N2:3.762',   # air components\n",
    "         'fuel':'CH3OCH3',          # fuel components\n",
    "         'type':'UV',               # sim type\n",
    "         'tot': 0.001,                # total time\n",
    "         }"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "ename": "CanteraError",
     "evalue": "\n***********************************************************************\nCanteraError thrown by ct2ctml_string:\nError converting input file \"mech/DME/DME_drg/DME_drg42.cti\" to CTML.\nPython command was: '/usr/local/opt/python/bin/python3.7'\nThe exit code was: 1\n-------------- start of converter log --------------\nTraceback (most recent call last):\n  File \"/usr/local/lib/python3.7/site-packages/cantera/ctml_writer.py\", line 2643, in convert\n    with open(filename, 'rU', **open_kw) as f:\nFileNotFoundError: [Errno 2] No such file or directory: 'mech/DME/DME_drg/DME_drg42.cti'\n\nDuring handling of the above exception, another exception occurred:\n\nTraceback (most recent call last):\n  File \"<stdin>\", line 9, in <module>\n  File \"/usr/local/lib/python3.7/site-packages/cantera/ctml_writer.py\", line 2670, in convert\n    text = open(filename, 'rU').readlines()\nFileNotFoundError: [Errno 2] No such file or directory: 'mech/DME/DME_drg/DME_drg42.cti'\n--------------- end of converter log ---------------\n***********************************************************************\n",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m-------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mCanteraError\u001b[0m                      Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-24-def147195dbe>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m      1\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mmech_i\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mmech_arr\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      2\u001b[0m     \u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmech\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mmech_i\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'name'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmech_i\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'path'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 3\u001b[0;31m     \u001b[0mgas\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mct\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mSolution\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mmech\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      4\u001b[0m     \u001b[0mdata\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mgetTcurv\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mgas\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mprops\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      5\u001b[0m     \u001b[0mplt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdata\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m't'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mdata\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'T'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32minterfaces/cython/cantera/base.pyx\u001b[0m in \u001b[0;36mcantera._cantera._SolutionBase.__cinit__\u001b[0;34m()\u001b[0m\n",
      "\u001b[0;32minterfaces/cython/cantera/base.pyx\u001b[0m in \u001b[0;36mcantera._cantera._SolutionBase._init_cti_xml\u001b[0;34m()\u001b[0m\n",
      "\u001b[0;31mCanteraError\u001b[0m: \n***********************************************************************\nCanteraError thrown by ct2ctml_string:\nError converting input file \"mech/DME/DME_drg/DME_drg42.cti\" to CTML.\nPython command was: '/usr/local/opt/python/bin/python3.7'\nThe exit code was: 1\n-------------- start of converter log --------------\nTraceback (most recent call last):\n  File \"/usr/local/lib/python3.7/site-packages/cantera/ctml_writer.py\", line 2643, in convert\n    with open(filename, 'rU', **open_kw) as f:\nFileNotFoundError: [Errno 2] No such file or directory: 'mech/DME/DME_drg/DME_drg42.cti'\n\nDuring handling of the above exception, another exception occurred:\n\nTraceback (most recent call last):\n  File \"<stdin>\", line 9, in <module>\n  File \"/usr/local/lib/python3.7/site-packages/cantera/ctml_writer.py\", line 2670, in convert\n    text = open(filename, 'rU').readlines()\nFileNotFoundError: [Errno 2] No such file or directory: 'mech/DME/DME_drg/DME_drg42.cti'\n--------------- end of converter log ---------------\n***********************************************************************\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAF+pJREFUeJzt3X+QXWWd5/H3t293pxMSSCAhYggT1LCKuIuaQSxnd1jZUaCmZCwdF3dnRIdaLFemZnatrcHZ3dLVoUprR62dGn8sDpRojQI7Y62pWbYodHDRmUEJgozAAD38MERMWhJISNM/7r3f/eOcTq7Y90dCd9/uc9+vqq4+9znPuXme9O3z6fM8z70nMhNJ0uAZ6ncDJEn9YQBI0oAyACRpQBkAkjSgDABJGlAGgCQNKANAkgaUASBJA8oAkKQBNdzvBnSycePG3LZtW7+bIUkryt133/2zzNzUrd6yDoBt27axa9eufjdDklaUiHiil3oOAUnSgDIAJGlAGQCSNKAMAEkaUAaAJA0oA0CSBpQBIEkDalm/D0DSi/P8TINv3LuHpw/PUG8kSdJMwFvBLnsvOWk1/+YNZyzqv2EASBX2xe88yqdve3jefRFL3Bgdk3O3rjcAJB2/x58+zCknjPI3V7+ZkdoQQwHhmV8lA0CqsMnpBqesHWVspNbvpmgZchJYqrAkCfyLX/PrGgARMRYR34+IH0bE/RHx38ryMyPiexExHhE3RcRoWb6qfDxe7t/W8lwfLssfioi3LlanJBUyHetXe71cAUwDb87MfwacC1wUEecDnwQ+k5mvAA4AV5T1rwAOlOWfKesREWcDlwGvBi4CPhcRXpdKi8i1PuqkawBk4bny4Uj5lcCbgb8oy28AfqPcvrR8TLn/wihmnS4FbszM6cx8DBgHzluQXkiaV3EF4CWA5tfTHEBE1CLiXmAfcBvwj8AzmVkvqzwJbCm3twC7Acr9zwKntJbPc4wkaYn1FACZ2cjMc4HTKf5qf+ViNSgiroyIXRGxa2JiYrH+GWlApFPAauuYVgFl5jPA7cAbgfURMbeM9HRgT7m9B9gKUO4/CXi6tXyeY1r/jWszc0dm7ti0qesdzSR14CSwOullFdCmiFhfbq8Gfg14kCII3llWuxz4Rrm9s3xMuf+vMzPL8svKVUJnAtuB7y9URyT9osQAUHu9vBHsNOCGcsXOEHBzZv5VRDwA3BgRfwTcA1xX1r8O+EpEjAP7KVb+kJn3R8TNwANAHfhgZjYWtjuSWmX6PgC11zUAMvM+4LXzlD/KPKt4MnMK+M02z3UNcM2xN1OStNB8J7BUYQ4BqRMDQKqwTBwAUlsGgFRhCV4CqC0DQKqwYhJYmp8BIEkDygCQKs4RILVjAEgV5iSwOjEApApL0k8DVVsGgFRhXgGoEwNAkgaUASBVmJ8Gqk4MAKnCvCm8OjEApArLxEkAtWUASBXm+V+dGACSNKAMAKnKnARWBwaAVGFOAqsTA0CqMJeBqhMDQKow7wimTgwASRpQBoBUYcUNYbwE0PwMAKnCHAJSJwaAVGGZ/W6BljMDQKqw4grASwDNzwCQpAFlAEhVlukUsNoyAKQKcxJYnRgAUoV5S0h1YgBIFeZN4dWJASBJA8oAkCrMISB1YgBIFeangaoTA0CqsOKNwCaA5mcASBWWmV4BqC0DQJIGlAEgVZwXAGrHAJAqzElgdWIASBXmTeHVSdcAiIitEXF7RDwQEfdHxO+V5R+NiD0RcW/5dUnLMR+OiPGIeCgi3tpSflFZNh4RVy9OlyTN8QpAnQz3UKcOfCgzfxAR64C7I+K2ct9nMvOPWytHxNnAZcCrgZcC34yIs8rdnwV+DXgSuCsidmbmAwvREUnSsekaAJn5FPBUuX0oIh4EtnQ45FLgxsycBh6LiHHgvHLfeGY+ChARN5Z1DQBpkfhpoOrkmOYAImIb8Frge2XRVRFxX0RcHxEbyrItwO6Ww54sy9qVS1pEzgGonZ4DICLWAn8J/H5mHgQ+D7wcOJfiCuFTC9GgiLgyInZFxK6JiYmFeEpJ0jx6CoCIGKE4+f95Zn4dIDP3ZmYjM5vAFzk6zLMH2Npy+OllWbvyn5OZ12bmjszcsWnTpmPtj6QW6V3h1UEvq4ACuA54MDM/3VJ+Wku1twM/Krd3ApdFxKqIOBPYDnwfuAvYHhFnRsQoxUTxzoXphqS2HAFSG72sAnoT8NvA30fEvWXZHwLvjohzKeaZHgfeD5CZ90fEzRSTu3Xgg5nZAIiIq4BbgRpwfWbev4B9kSQdg15WAX2X+f+GuKXDMdcA18xTfkun4yQtLAeA1InvBJYqzhEgtWMASNKAMgCkKnMMSB0YAFLFhW8FVhsGgCQNKANAqjBHgNSJASBVnANAascAkCrMj4JQJwaAVHHOAasdA0CSBpQBIFWYA0DqxACQKs4RILVjAEjSgDIApApzEZA6MQCkivOjINSOASBJA8oAkCosXQekDgwAqeIcAFI7BoAkDSgDQKowVwGpEwNAqjrHgNSGASBJA8oAkCrMISB1YgBIFReOAakNA0CSBpQBIEkDygCQKs6PAlI7BoAkDSgDQKowbwqvTgwAqeIcAVI7BoAkDSgDQKowB4DUiQEgVZyrgNSOASBJA8oAkCrMRUDqxACQKs7PAlI7BoAkDSgDQKowbwqvTgwAqeJcBaR2ugZARGyNiNsj4oGIuD8ifq8sPzkibouIR8rvG8ryiIg/iYjxiLgvIl7X8lyXl/UfiYjLF69bkqRuerkCqAMfysyzgfOBD0bE2cDVwLcyczvwrfIxwMXA9vLrSuDzUAQG8BHgDcB5wEfmQkPS4nAVkDrpGgCZ+VRm/qDcPgQ8CGwBLgVuKKvdAPxGuX0p8OUs3Amsj4jTgLcCt2Xm/sw8ANwGXLSgvZH0CxwCUjvHNAcQEduA1wLfAzZn5lPlrp8Cm8vtLcDulsOeLMvalUuS+qDnAIiItcBfAr+fmQdb92XxmbMLcrEZEVdGxK6I2DUxMbEQTykNLEeA1ElPARARIxQn/z/PzK+XxXvLoR3K7/vK8j3A1pbDTy/L2pX/nMy8NjN3ZOaOTZs2HUtfJL1AZjLkGJDa6GUVUADXAQ9m5qdbdu0E5lbyXA58o6X8PeVqoPOBZ8uholuBt0TEhnLy9y1lmaRF0kwMALU13EOdNwG/Dfx9RNxblv0h8Ang5oi4AngCeFe57xbgEmAcmATeB5CZ+yPi48BdZb2PZeb+BemFpHk1Mxny/K82ugZAZn6X9jcVunCe+gl8sM1zXQ9cfywNlHT8ms0kvAJQG74TWKqwdAhIHRgAUoU5BKRODACpwpoJQyaA2jAApAprZvpOYLVlAEgV5hyAOjEApApzDkCdGABShTV9J7A6MACkCmsmvg9AbRkAUkVleTMAh4DUjgEgVVSz/ChQh4DUjgEgVVTTKwB1YQBIFTUXAM4BqB0DQKqodAhIXRgAUkU5BKRuDACpopwEVjcGgFRRR+cA+twQLVsGgFRR2Sy+ewWgdgwAqaK8AlA3BoBUUY0yAIadBVYbBoBUUfVGEQC1IX/NNT9fGVJF1ZvFJMBwzSsAzc8AkCpq7grAISC1YwBIFVUv3wgwXPPXXPPzlSFV1JEhIK8A1IYBIFWUQ0DqxgCQKqpxZAjIAND8DACpoo4OAflrrvn5ypAqyiEgdWMASBXlKiB14ytDqqi5AKh5BaA2DACpouqNYg5gxElgtWEASBXlFYC6MQCkijqyDNRVQGrDV4ZUUbMNPwxOnRkAUkUdvQIwADQ/A0CqqCPvA3AZqNrwlSFVVN0rAHVhAEgV5aeBqhsDQKqoox8F4a+55tf1lRER10fEvoj4UUvZRyNiT0TcW35d0rLvwxExHhEPRcRbW8ovKsvGI+Lqhe+KpFbeElLd9PKnwZeAi+Yp/0xmnlt+3QIQEWcDlwGvLo/5XETUIqIGfBa4GDgbeHdZV9Ii8Y1g6ma4W4XMvCMitvX4fJcCN2bmNPBYRIwD55X7xjPzUYCIuLGs+8Axt1hST/w0UHXzYgYHr4qI+8ohog1l2RZgd0udJ8uyduW/ICKujIhdEbFrYmLiRTRPGmwz9Sa1oXAZqNo63lfG54GXA+cCTwGfWqgGZea1mbkjM3ds2rRpoZ5WGjgzjSajnvzVQdchoPlk5t657Yj4IvBX5cM9wNaWqqeXZXQol7QIZupNRocNALV3XK+OiDit5eHbgbkVQjuByyJiVUScCWwHvg/cBWyPiDMjYpRionjn8TdbUjfTBoC66HoFEBFfAy4ANkbEk8BHgAsi4lwggceB9wNk5v0RcTPF5G4d+GBmNsrnuQq4FagB12fm/QveG0lHTNcbrDIA1EEvq4DePU/xdR3qXwNcM0/5LcAtx9Q6ScfNISB146tDqqiZupPA6sxXh1RR0/WmQ0DqyFeHVFEz9Sarhmv9boaWMQNAqqiZhnMA6sxXh1RRz880GBvxV1zt+eqQKurwTJ21q47rvZ4aEAaAVFGHp+ucYACoAwNAqqhDU14BqDMDQKqg2UaT6XrTAFBHBoBUQYen6wAOAakjA0CqoENTRQCsHTMA1J4BIFXQxHPTAGxau6rPLdFyZgBIFfSzQ2UArDMA1J4BIFXQ3BXARq8A1IEBIFXQvoNFAJyydrTPLdFyZgBIFfTj/ZNsWb+aET8OWh346pAq6NGfHWbbxjX9boaWOQNAqpjM5LGJ59h2ygn9boqWOQNAqpgnnp7k4FSds196Yr+bomXOAJAq5p7dBwB43Rkb+twSLXcGgFQxfzv+NOvGhjlr87p+N0XLnAEgVUi90eSbD+7lwleeSm0o+t0cLXMGgFQh335oggOTs1x0zmn9bopWAANAqpA/++6jnHbSGBe+6tR+N0UrgAEgVcTt/7CPOx/dzxW/cqZvAFNPfJVIFfDs5Cz/9Rs/4mWbTuA9b9zW7+ZohfDDwqUVbrre4HdvvIe9B6e46f1vZHTYv+vUGwNAWsEmZ+pc9dV7uOPhCT75jte49l/HxACQVqhH9h7iqq/ew8P7DnHN28/hX//yGf1uklYYA0BaYQ5NzfKnt49z/XcfY93YCF/+nfP459s39btZWoEMAGmF2Hdwii/97eN85c4nODRV5zdffzp/cPErvemLjpsBIC1jh6frfPPBvXz9B3v4ziMTJHDxOS/hA7/6Cl5z+kn9bp5WOANAWkYazeQffnqQ7zzyM/7fQxPsemI/s41ky/rVfOCCl/PO12/lzI1+zLMWhgEg9Ulm8pNnp3jgJwe5d/cB7vnxM/xw9zMcnmkA8MqXrON3fuVM3vxPTuWXt53MkJ/towVmAEiLbGq2wZ5nnufH+yd5dOIwD//0EA/vO8Qje5/juek6AMNDwatOO5F3vP50zt26nje9YiObTxzrc8tVdQaA9CI0m8mByRn2HZpm78Ep9h2aZs+B59l9YJLd+yfZvf95fnpw6ueOOfmEUc7avJZ3vG4L2zev45UvWcerX3oSq0drfeqFBpUBILXITJ6fbXBgcpYDh2c4MDnDgclZnpmcYf/hGZ5+boZ9h6bYe3CafQenmHhumtlG/txzRMBLThxj68lreNMrNnLGyWvYevJqzjh5Dds2nuCqHS0bBoAqZe4E/txUnYNTdZ6brnNoapZDU/WyrNxuKX9mcpYDkzM8MznL/skZZurNts+/fs0Ip65bxeYTx3jZplPYfOLYkcdHvp+4ilXD/jWv5a9rAETE9cCvA/sy85yy7GTgJmAb8Djwrsw8EBEB/A/gEmASeG9m/qA85nLgv5RP+0eZecPCdkUrSb3R5PnZBs/PNJicafD8bPl9psHkTP3I46Ks3rLdYHL2F8sOzxQn9eem6tSb2fXfP2G0xrqxEdaODbN+9QhbT17DPz19hA1rRtlwwigb1oywfs0oJ7dsr189wrCfsqkK6eUK4EvAnwJfbim7GvhWZn4iIq4uH/8BcDGwvfx6A/B54A1lYHwE2AEkcHdE7MzMAwvVEfWm2UxmGk1mG01mG8lso8lMvXm0rF7sn643mK43mZ5t3S6/v3C73mBqrt7s0bK546eOlB89ptHDSbrV8FCwerTGmtEaq0dqrB4dZs1ojXVjw5y6bhUnrBpm3dgwa1cNs25shHVjwy1fI2X50W3vliX1EACZeUdEbHtB8aXABeX2DcC3KQLgUuDLmZnAnRGxPiJOK+velpn7ASLiNuAi4Gsvugd90mgm9WaTeiOLr2aTerM4oTaayWyjZX8zqTea5fdkttmkUe6fqzd34p1tOTkfOTHXi7KZ8oQ9d9I+UlY/etxMI5mpN46c3I/WLZ7/WE+8nYyNDLFquMaq4SFWtW4PDzE2UmPtquGibGSoLK8d2Td3zOrRYdaMlCf20RpryhP73Ml+zcgwq0drfsKltAiOdw5gc2Y+VW7/FNhcbm8BdrfUe7Isa1e+KJ6ZnOGdX/g7MpMEMoux4WZCkuVjjuxvZlFWnBvnto8e28yEhHozi5N7s0ku3Hm0o5FaMFIbYnR4qPheGzpSNlIbYmR4iFW14mS6dmz45+rMHXP0+KPPNVo7uq+17mhtiOFasGq4dvQE33oCL7dHa0MUI36SVqoXPQmcmRkRC3Y6jIgrgSsBzjjj+D7dsDYUnLV5LUEQARHBUEAUz1+UlfuGyu2hIeAFZcV2lO0qhiFqQ8UJc3ioOFEODwXDtaHyezAyNERtbrtWbB+pP1d37rjyuWpDceSk3HrSHqmFJ1lJi+Z4A2BvRJyWmU+VQzz7yvI9wNaWeqeXZXs4OmQ0V/7t+Z44M68FrgXYsWPHcQXLurERPvdvX388h0rSwDjegdWdwOXl9uXAN1rK3xOF84Fny6GiW4G3RMSGiNgAvKUskyT1SS/LQL9G8df7xoh4kmI1zyeAmyPiCuAJ4F1l9VsoloCOUywDfR9AZu6PiI8Dd5X1PjY3ISxJ6o/IpZrNPA47duzIXbt29bsZkrSiRMTdmbmjWz3X1knSgDIAJGlAGQCSNKAMAEkaUAaAJA2oZb0KKCImKJaZHq+NwM8WqDnLmf2slkHpJwxOX5e6n7+UmZu6VVrWAfBiRcSuXpZCrXT2s1oGpZ8wOH1drv10CEiSBpQBIEkDquoBcG2/G7BE7Ge1DEo/YXD6uiz7Wek5AElSe1W/ApAktbHiAyAiLoqIhyJivLw/8Qv3r4qIm8r935vn9pYrRg99/Y8R8UBE3BcR34qIX+pHO1+sbv1sqfeOiMiIWHarK3rRSz8j4l3lz/T+iPjqUrdxIfTwuj0jIm6PiHvK1+4l/WjnixUR10fEvoj4UZv9ERF/Uv4/3BcRr1vqNv6CzFyxX0AN+EfgZcAo8EPg7BfU+ffAF8rty4Cb+t3uRezrvwTWlNsfWIl97aWfZb11wB3AncCOfrd7kX6e24F7gA3l41P73e5F6ue1wAfK7bOBx/vd7uPs678AXgf8qM3+S4D/S3FzwvOB7/W7zSv9CuA8YDwzH83MGeBGihvTt7qU4sb1AH8BXBgr8z6LXfuambdn5mT58E6KO6+tNL38TAE+DnwSmFrKxi2gXvr574DPZuYBgMzcx8rTSz8TOLHcPgn4yRK2b8Fk5h1Ap/ucXAp8OQt3AuvLOyr2zUoPgF5uNn+kTmbWgWeBU5akdQurl762uoLir42Vpms/y0vnrZn5f5ayYQusl5/nWcBZEfE3EXFnRFy0ZK1bOL3086PAb5U3nLoF+N2ladqSO9bf4UX3om8Kr+UnIn4L2AH8ar/bstAiYgj4NPDePjdlKQxTDANdQHE1d0dEvCYzn+lrqxbeu4EvZeanIuKNwFci4pzMbPa7YVW30q8A2t2Eft46ETFMcYn59JK0bmH10lci4l8B/xl4W2ZOL1HbFlK3fq4DzgG+HRGPU4yl7lyBE8G9/DyfBHZm5mxmPgY8TBEIK0kv/bwCuBkgM/8OGKP47Jyq6el3eCmt9AC4C9geEWdGxCjFJO/OF9RpvYH9O4G/znJGZoXp2teIeC3wPylO/itxvBi69DMzn83MjZm5LTO3Ucx1vC0zV9q9Q3t57f5vir/+iYiNFENCjy5lIxdAL/38MXAhQES8iiIAJpa0lUtjJ/CecjXQ+cCzmflUPxu0ooeAMrMeEVcBt1KsNrg+M++PiI8BuzJzJ3AdxSXlOMUEzWX9a/Hx67Gv/x1YC/yvcp77x5n5tr41+jj02M8Vr8d+3gq8JSIeABrAf8rMFXX12mM/PwR8MSL+A8WE8HtX4h9pEfE1isDeWM5nfAQYAcjML1DMb1wCjAOTwPv609KjfCewJA2olT4EJEk6TgaAJA0oA0CSBpQBIEkDygCQpAFlAEjSgDIAJGlAGQCSNKD+P/55LShFj6GqAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "for mech_i in mech_arr:\n",
    "    name, mech = mech_i['name'], mech_i['path']\n",
    "    gas = ct.Solution(mech)\n",
    "    data = getTcurv(gas, props)\n",
    "    plt.plot(data['t'],data['T'])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
